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Obligate symbiotic fungi that form arbuscular mycorrhizae (AMF; belonging to the 
Glomeromycota phylum) are some of the most important soil microorganisms. AMFs 
facilitate mineral nutrient uptake from the soil, in exchange for plant-assimilated carbon, 
and promote water-stress tolerance and resistance to certain diseases. AMFs colonize 
the root by producing inter- and intra-cellular hyphae. When the fungus penetrates the 
inner cortical cells, it produces a complex ramified structure called arbuscule, which is 
considered the preferential site for nutrient exchange. Direct DNA extraction from the 
whole root and sequencing of ribosomal gene regions are commonly carried out to 
investigate intraradical AMF communities. Nevertheless, this protocol cannot discriminate 
between the AMFs that actively produce arbuscules and those that do not. To solve 
this issue, the authors have characterized the AMF community of arbusculated cells 
(AC) through a laser microdissection (LMD) approach, combined with sequencing-based 
taxa identification. The results were then compared with the AMF community that was 
found from whole root DNA extraction. The AMF communities originating from the LMD 
samples and the whole root samples differed remarkably. Five taxa were involved in the 
production of arbuscules, while two taxa were retrieved inside the root but not in the 
AC. Unexpectedly, one taxon was found in the AC, but its detection was not possible 
when extracting from the whole root. Thus, the LMD technique can be considered a 
powerful tool to obtain more precise knowledge on the symbiotically active intraradical 
AMF community. 
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INTRODUCTION 

Obligate symbiotic fungi that form arbuscular mycorrhizae 
(AMF; belonging to the phylum Glomeromycota) are among the 
most important soil microorganisms. AMFs facilitate mineral 
nutrient uptake from the soil in exchange for plant-assimilated 
carbon and promote water-stress tolerance and resistance to cer- 
tain diseases (Smith and Read, 2008). About 240 AMF species 
have been described so far (Oehl et al., 2011; Kriiger et al., 2012), 
primarily on the basis of spore morphology but increasingly 
supported by DNA sequence information. The use of next gen- 
eration sequencing methods (e.g., pyrosequencing) has changed 
our understanding of the diversity of AMF, by allowing AM fungi 
to be directly detected from environmental samples (Opik et al., 
2009, 2013; Lumini et al, 2010; Dumbrell et al, 2011; Davison 
et al., 2012). Comparisons between AMF communities inside and 
outside the roots have shown that all the phylotypes from the 
roots were represented in the soil, which, as expected, maintains 
a higher AMF biodiversity than roots (Balestrini et al., 2010). 

The AMF enters the root epidermis through the hyphopodium 
and starts to colonize the root by producing inter and intracellu- 
lar hyphae. When the fungus penetrates the inner cortical cells, 
it produces a complex ramified structure called arbuscule, which 
is considered the preferential site for nutrient exchange (Smith 



and Read, 2008). It is rather easy, although quite time-demanding 
and subjective, to spot arbuscules inside root parenchyma cortex 
cells through staining and microscopy (Trouvelot et al., 1986). 
Unfortunately, it is rather difficult, if not impossible, to dis- 
criminate the specific taxa involved in the production of these 
structures according to the morphology, even at a high tax- 
onomical rank level. In order to investigate intraradical AMF 
communities, direct DNA extraction from the whole root and 
sequencing of ribosomal gene regions are commonly carried out 
(Balestrini et al., 2010; Lumini et al., 2011; Borriello et al., 2012; 
Pellegrino et al., 2012). Nevertheless, this protocol is not able 
to discriminate between the AMFs that actively produce arbus- 
cules and those that do not. Apart from arbuscules, a multitude 
of intraradical AMF structures are normally encountered in a 
host-root, i.e., hyphae, vesicles, spores, and hyphopodia, and they 
contain DNA but have no direct correlation with the plant-fungus 
exchange activity. Therefore, the amplification of DNA from such 
structures alters the actual shape of the active AMF community 
due to a dilution effect. 

In an attempt to discriminate the active intraradical AMF 
community, the authors have proposed the characterization 
of the AMF community of arbusculated cells (AC) through 
a laser microdissection (LMD) approach combined with 
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sequencing-based taxa identification. The idea of applying LMD 
to plant tissues was first advanced during the early years of the 
new millennium and its potential was reviewed by Day et al. 
(2005). As predicted by these authors, the LMD technique has 
resulted to be a powerful tool for isolating specific tissues and cell 
types from sectioned plant specimens, as it allows the cell spe- 
cific extraction of RNA, DNA, and proteins/metabolites (Nelson 
et al., 2006; Balestrini and Bonfante, 2008). During the last few 
years, LMD has been used to study gene expression in AMF. 
Particular attention has been paid to specific plant and fun- 
gal genes expressed in the cortical cells containing arbuscules 
(Balestrini et al, 2007; Gomez et al, 2009; Guether et al., 2009; 
Gaude et al., 2011; Hogekamp et al., 2011; Perez-Tienda et al., 
2011; Tisserant et al., 2012). However, no information is available 
on the potential use of LMD for the identification of intraradical 
active AMF taxa on the basis of only DNA extraction from homo- 
geneous cell type populations (AC). In the present work, the aim 
has been to isolate the specific fungal compartments that have a 
crucial symbiotic role in order to correctly identify functionally 
active AMF species. 

In this study, Camellia japonica L. 'Nuccio's Pearl' was cho- 
sen as the experimental plant. This plant species is a broadleaf 
evergreen perennial shrub and it is difficult to consider it a 
model plant or agricultural crop, since very limited research 
information is available on the species. The camellia is an aci- 
dophilic plant and is characterized by a very specific rhizo- 
spheric chemistry because of the phenolic-rich root exudates 
and low pH characteristics that make this species a difficult 
plant material to work with. The rhizosphere of C. japonica 
and its close species has been demonstrated to host a diver- 
sity of microorganisms, including AMFs (Gupta and Sharma, 
2010), bacteria (Pandey and Palni, 1996) and nematodes (Zheng 
et al., 2000). The authors have developed a method to extract 
DNA from microdissected arbuscule-containing cells from the 
mycorrhizal roots of potted camellias. After DNA amplifica- 
tion with AMF specific primers, the PCR products were cloned 
and sequenced. This led to the identification of the AM fungi 
that formed arbuscules inside the roots. These data were com- 
pared with the DNA sequences obtained from direct whole 
root DNA extraction in order to provide information on the 
AMF species that are involved in the formation of a functional 
symbiosis. 

MATERIALS AND METHODS 

PLANT MATERIAL AND GROWTH CONDITIONS 

C. japonica 'Nuccio's Pearl' was selected as the test plant because 
of its commercial and ornamental value. The cultivation cycle 
lasted almost 2 years (February 2011-November 2012) and was 
carried out in a commercial nursery devoted to the production 
of acidophilic ornamental plants (Tecnoverde s.p.a., Verbania, 
Piedmont — Italy). A substrate mixture with 50% of coconut 
fiber and 50% of standard substrate (89% commercial Sphagnum 
peat and 11% agriperlite) was used as the growing medium for 
all the plants. This substrate had a pH of 5.6 (substrate/water 
ratio: 1/5 v/v), an electrical conductivity (EC) of 38 u,S/cm and 
a C/N ratio of 76.7. The K, Ca, and Mg contents were 0.1, 
0.4 and 0.1%, respectively, while no detectable concentrations 



of P were present. Neither a spore count nor a PCR-based 
approach with specific primers (see the Materials and Methods 
section below) showed the presence of AMF propagules in the 
substrate mixture. A fertilizer dosage of 2.25 kg m~ 3 CaC03, 
1.12 kg nr 3 Osmocote Exact® (15N-3.9P-7.5K + trace micronu- 
trients, 8-9 months, Scotts, Merysville, Ohio) and 0.5 kg m~ 3 
Scorie Thomas (6.9P, phosphatic inorganic fertilizer, Timac Italia, 
Ripalta Arpina, Italy) were added to the substrate mixture. In 
the first week of February 2011, three camellia rooted cuttings 
were planted per pot (9 cm in diameter) for a total of 35 pots 
(105 rooted cuttings), and 100 g of AMF inoculum was added 
to each pot (see below for the methods regarding the inocu- 
lum formulation). The pots were randomly distributed over a 
bench in a frost-free greenhouse (2.5°C minimum, 16°C aver- 
age). The shading was set to a constant 30%. 40% additional 
shading was provided when the natural illumination exceeded 
60 Klux. No supplementary lighting was applied. All the camel- 
lias were watered and fertigated with 20N-8.7P-16.6K Peter 
Professional hydrosoluble fertilizer (Scotts, Merysville, Ohio) at 
0.8-0.9 g 1~' about once every 2 weeks from March to October, 
according to the weather conditions. After 9 months of culti- 
vation, the plants were transplanted into 15 cm diameter pots 
and pruned. The cultivation ended in November 2012, when the 
plant roots were harvested and stored at — 80° C or immediately 
processed. 

INOCULUM PREPARATION 

The bulk soil surrounding centennial camellia specimens in 
the public garden on the Isola Madre isle (45°54'43.46" N, 
8°32'21.35" E, Verbania, Italy) was collected on May 22nd 2009. 
A physicochemical analysis was performed in three replicates to 
characterize the soil properties, according to European Standard 
methods (EN). The soil collected had a pH of 5.9, an EC of 
62.7 (iS/cm, and very low porosity (9%). The C/N ratio was 13.1, 
the total P was 0.02% while the K, Ca, and Mg were 0.45, 0.16, 
and 4.86%, respectively. The sampling operations consisted in 
digging to the first 5-20 cm and collecting bulk soil (ca 1.5 Kg) 
and fine feeder roots belonging to a centennial C. japonica Alba 
Plena' tree. The bulk soil and roots were taken to the lab, gently 
homogenized and mixed with 60% sterile sand. The resulting 
substrate was used as growing media for trap plants in order 
to propagate the AMF portion present in the inoculum. The 
species used for this purpose were Nicotiana tabacum, Lolium 
perenne, Trifolium repens, and Plantago major, all of which are 
recognized as being highly efficient trap plants (Liu and Wang, 
2003; Oehl et al, 2004). After 9 months, the AMF inoculum 
was re-propagated once again by homogenizing the trap plant 
substrate and mixing it with 60% sterile sand. The amount of 
inoculum required to carry out the experiment was obtained after 
another 9 months. The total inoculum, including chopped root 
systems, was homogenized and air-dried at room temperature. 
A portion of roots belonging to the trap plant species was 
washed free of the substrate, air-dried at room temperature and 
immediately used for a morphological assessment of mycorrhiza- 
tion levels, according to the protocol described by Trouvelot 
et al. (1986). The mycotrophic status of the root samples was 
consistent, the frequency of mycorrhization was more than 85%, 
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the intensity of colonization was around 30%, and the presence 
of arbuscules in the mycorrhizal roots amounted to 70%. In order 
to evaluate the infectiveness of the inoculum, the Most Probable 
Number (MPN, Porter, 1979) of AMF propagules was calculated 
at 130 propagules/cm 3 . 

MOLECULAR CHARACTERIZATION OF THE AMF INOCULUM 

Before starting the experimentation on the camellia plants, a 
molecular AMF community composition study was performed 
on the crude AMF inoculum, in order to characterize all the 
isolates that had propagated successfully. DNA extraction was 
performed on the trap plant roots and on the crude inoculum, 
as reported by Borriello et al. (2012), according to the proto- 
cols for frozen root samples and bulk soil, respectively. A nested 
PCR approach (Borriello et al., 2012) was adopted to amplify 
a region of the small subunit (SSU) of the Glomeromycota 
ribosomal DNA. This approach consisted of a first amplifi- 
cation with the universal eukaryotic primers NS1 and NS4 
(White et al., 1990) and a following amplification round on 
the 100X diluted PCR product with Glomeromycota-specific 
primers AML1 and AML2 (Lee et al, 2008). PCR was car- 
ried out using 0.2 mM dNTPs, 3.5 mM of MgCl 2 , 0.5 [iM of 
each primer and the supplied reaction buffer, with two units 
of High Fidelity Taq (ROCHE) to obtain a final volume of 
20 (jlI. Amplifications were carried out in 0.2 ml PCR tubes 
using a Biometra T Gradient thermocycler, according to the fol- 
lowing steps: 5min initial denaturation at 95° C; 35 cycles of 
lmin at 95°C, 1 min at 55°C and 58°C for the two nested 
PCR rounds, respectively, 1 min at 72° C; and a final elonga- 
tion of 10 min at 72° C. A negative control was included in 
the PCR to check for contamination. All the PCR products 
were checked using 1.2% agarose gel stained with Ethidium 
Bromide (Sigma-Aldrich). The PCR products were purified using 
QIAquick (Qiagen, Hilden, Germany), cloned in a pGEM-T 
Easy vector (Promega, Madison, Wisconsin, USA) and trans- 
formed into competent Escherichia coli cells (Xll blue). After 
colony PCR, inserts of approximately the correct expected size 
(~800bp) were submitted to restriction fragment length poly- 
morphism (RFLP) analysis using two endonucleases, Hsp92II 
(|CATG|) and Hinfl (G|ANT|C), in order to reduce the num- 
ber of clones to be sequenced. Subsequently, 1-2 clones were 
selected as being representative of each found RFLP profile and 
were sequenced using T7 vector primers by Macrogen sequencing 
services (Macrogen, Seoul, Korea). Sequence editing was con- 
ducted using SEQUENCER V4.2.2 (Gene Codes Corporation, 
Ann Arbor, MI, USA). The potential chimera sequences were 
identified using the Chimera Detection programme (Cole et al., 
2003). A search for similar sequences was conducted with the 
BLAST tool (Zhang et al, 2000), provided by GenBank. All the 
sequences were aligned using the multiple sequence comparison 
alignment tool by MAFFT version 6 (Katoh and Toh, 2008) and 
a Neighbor-joining (NJ) phylogenetic analysis (5000 bootstrap 
replicates) was performed with MEGA version 5 (Tamura et al., 
201 1). Sequences of the Paraglomerales AMF order were set as the 
outgroup. A Maximum Likelihood phylogenetic analysis (1000 
bootstrap replicates) was also conducted and comparable results 
were obtained. 



ASSESSMENT OF THE INOCULATED ROOT COLONIZATION BY AM 
FUNGI 

A portion of inoculated camellia roots from three randomly 
selected potted plants was washed free of the substrate, air-dried 
at room temperature and chopped into small fragments that were 
cleared twice with 10% KOH for 30 min at 90° C. After rinsing 
several times with tap water, the root fragments were trans- 
ferred to alkaline water (H20:NH 4 OH:H 2 02 198:1:1) for 60 min, 
stained with 0.1% cotton blue in 80% lactic acid for about 18 h 
and then de-stained several times with 80% lactic acid washes. 
The roots were cut into 1 cm fragments and placed onto micro- 
scope slides. 25 fragments were observed for each replicate, for a 
total of 75 root fragments. Mycorrhiza frequency (F%), root AMF 
colonization intensity in the mycorrhizal root parts (m%) and in 
the whole root system (M%), and the presence of arbuscules in 
mycorrhizal root parts (a%) and in the whole root system (A%) 
were determined and calculated, as described by Trouvelot et al. 
(1986). 

DNA EXTRACTION FROM THE WHOLE ROOT AND LASER 
MICRODISSECTED CELLS 

DNA extraction on the whole root was performed in duplicate 
(abbreviated as sample Rl and R2) with the DNeasy Plant Mini 
Kit (Qiagen, Crawley, UK), on 0.1 g of finely ground frozen root 
samples taken from three randomly chosen camellia pots. In order 
to prepare the plant material for microdissection, the roots from 
the same camellia pots were fixed in freshly prepared absolute 
ethanokacetic acid (3:1) at 4°C overnight for paraffin embed- 
ding (Balestrini et al, 2007). A Leica AS LMD system (Leica 
Microsystems) was used to isolate the arbuscule-containing cells 
from the prepared tissue sections, as described in previous works 
(Balestrini et al, 2007; Guether et al, 2009). Three different 
aliquots containing 100, 300, and 500 putative AC (abbreviated 
as A100, A300, and A500, rispectively) were collected in 0.2 ml 
eppendorf tube caps from slides filled with root transversal and 
longitudinal ll-12|xm root sections obtained from at least 10 
different root fragments (Figure 1). After collection, the sam- 
ples were briefly centrifuged to move the plant material at the 
bottom of the tubes and a DNA extraction buffer (including 
Proteinase K) from the Arcturus PicoPure DNA Extraction kit 
(Life Technologies) was added. The samples were gently vor- 
texed and incubated at 65° C for 3h. In order to inactivate 
Proteinase K, the samples were placed at 95°C for 10 min. After 
DNA extraction, yield was measured using a NanoDrop 1000 
spectrophotometer. 

NESTED PCR. CLONING, RFLP. AND SEQUENCING OF THE FUNGAL 
RIBOSOMAL (rRNA) GENE 

Nested PCR, cloning, and RFLP typing were carried out for whole 
root or microdissected arbusculated cell samples with the same 
approach adopted for the inoculum molecular characterization. 
Two microliter of the undiluted DNA was used as input during the 
first round of nested PCR for the microdissected arbusculated cell 
samples and 2 \il of the 1:20 diluted DNA was used for the whole 
root samples. The only modification consisted of the use of an 
undiluted PCR product as a template for the second round of PCR 
when the microdissected samples were considered. The same way 
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FIGURE 1 | Sections of different camellia roots on a Leica RNase-free 
PEN-foil slide after paraffin embedding. Arbuscule-like structures located 
on a transversal root section (A). Transversal root section before (B) and after 



(C) multiple contiguous arbuscules were collected using LMD. Transversal 
root section before (D,E) and after (F) a single arbuscule was collected 
using LMD. 



as during the characterization of the inoculum, after colony PCR, 
inserts of approximately the correct size from each SSU DNA 
library were RFLP typed by comparing the obtained bandage 
with the RFLP profile data found during the previous molecu- 
lar characterization of the inoculum. 24 representative samples 
were sequenced, using primer pair AML1-AML2, by Sequencing 
Service at Munich University (Database ID: http://www.gi.bio. 
lmu.de/sequencing). The sequences were edited and their affili- 
ation investigated, as described in the above sections. In order to 
assess the reliability of the community analysis, the representa- 
tiveness of the observed community distribution was calculated, 
with respect to the real community distribution in the sample, 
with the Estimated Sample Coverage (ESC) and the bias corrected 
CHAOl species richness estimator, by using the freeware SPADE 
(Chao and Shen, 2010). In addition, rarefaction curves were built 
with PAST version 2.17 (Hammer et al, 2001). 

RESULTS 

AMF INOCULUM COMPOSITION 

Fifty-three transformed clones were sequenced as representatives 
of each of the 35 RFLP profiles that were found. The phylogenetic 
inference (Figure 2) revealed that, of the 35 obtained RFLP pro- 
files, 27 (45 sequences) clustered with major taxonomical groups 
belonging to the Glomeromycota phylum, i.e., Archaeosporales, 
Glomerales, and Diversisporales. The remaining 8 RFLP pro- 
files (8 sequences) were mainly related to sequences from plant 
and other soil organisms, as they were the result of aspecific 
amplifications. 

MYC0TR0PHIC STATUS OF INOCULATED CAMELLIA ROOTS 

Most of the camellia root fragments exhibited a very high degree 
of AMF colonization. The mean percentage value of the frequency 
was 90.0 ± 7.3, while the percentage intensity of the coloniza- 
tion calculated in the root system and in the root fragments 
was 57.1 ± 6.2 and 63.1 ± 1.7, respectively. Arbuscules covered 



54.9 ± 5.6% of the total root system. In particular, 96.4 ± 2.2% of 
the mycorrhizal parts hosted these structures. Typical symbiotic 
scenarios are portrayed in Figure 3. 

AMF DETECTION IN THE ARBUSCULATED CELLS AND WHOLE ROOT 
SAMPLES 

The DNA extractions yielded 25.8, 29.2, and 27.9 ng ijlI 1 for 
the LMD samples (A100, A300, A500) and 20.5, and 26.0 ng 
jxl 1 for the whole root samples (Rl and R2), respectively. 
The nested PCR, using the AML1-AML2 primer combination, 
yielded a clone library that contained SSU fragments from 
the whole root and arbusculated cell matrices for each sample 
under study. 120 clones were screened and RFLP typed. The 
presence of 13 different RFLP profiles was highlighted and 24 
representative samples (1-3 per RFLP profile) were sequenced. 
Apart from sequences showing high similarity (>97% iden- 
tity) with AM fungi (Glomeromycota phylum), sequences corre- 
sponding to non-target organisms were also obtained, as already 
shown during the inoculum molecular characterization. BLAST 
searches revealed that, of the 24 obtained SSU sequences, 16 
showed high similarity to sequences from taxa belonging to the 
Glomeromycota phylum, while the remaining eight were mainly 
related to sequences from other organisms (e.g., C. japonica and 
other fungi) and were the result of an aspecific amplification of 
primer pair AML 1 - AML2 . The percentage of screened clones con- 
taining a non-AMF SSU DNA fragment is shown in Figure 4. 
None of the clones screened from the SSU DNA library belonging 
to the aliquot containing 100 arbusculated cells (A100) contained 
an AMF DNA fragment. The DNA extracted from the aliquots 
containing a higher number of cells, 300 (A300) and 500 (A500), 
showed a proportional increase in terms of AMF SSU fragments 
(around 50 and 70%, respectively). The SSU DNA libraries from 
the whole root samples (Rl and R2) showed little or no pres- 
ence of non-AMF fragments. After discarding all the non-target 
fragments, clones showing the same RFLP profiles were affiliated 
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FIGURE 2 | The phylogenetic tree constructed on a portion of the 
arbuscular mycorrhizal fungal SSU rDNA (~800 bp) using the 
Neighbour-Joining method (Saitou and Nei, 1987). The percentage of 
replicate trees in which the associated taxa clustered together in the 
bootstrap test (5000 replicates) is shown next to the branches (Felsenstein, 
1985). Branches corresponding to partitions reproduced in less than 45% 
bootstrap replicates were collapsed. The evolutionary distances were 



computed using the Kimura 2-parameter method (Kimura, 1980) and are in 
the units of the number of base substitutions per site. The analysis involved 
79 nucleotide sequences (52 reference sequences and 27 sequences 
representative of the found RFLP profiles, highlighted in bold black). 
Evolutionary analyses were conducted in MEGA5 (Tamura et al., 2011). The 
external brackets indicate diverse taxonomic groups, as described by Kruger 
et al. (2012). The Paraglomerales order was used as the outgroup. 



to an AMF phylogenetic group on the basis of BLAST search 
results, following the same procedure adopted during the inocu- 
lum molecular characterization. As a result, an AMF community 
distribution was obtained for each sample (Figure 5). The AMF 
community of sample A300 was characterized by 12 observed 
individuals and the presence of 3 RFLP profiles, RHIZO-1, the 
most abundant (66.7%), GLO-6 (25%), and RHIZO-2 (8.3%). 
The ESC for this sample was 91.7% but the rarefaction curve 
did not reach the asymptote (Figure 6), while bias corrected 
CHAOl (3.0 ±0.4) confirmed the reliability of the number 
of found RFLP profiles. In A500, the community counted 23 
observed individuals distributed over 4 RFLP profiles, ARCH-1, 



the most abundant (52.2%), ARCH-2 (30.4%), RHIZO-1 (8.7%), 
and RHIZO-2 (8.7%). The sample coverage was maximum (ESC 
= 100%), the related rarefaction curve reached the asymptote 
(Figure 6) and bias corrected CHAOl (4.0 ± 0.0) confirmed the 
representativeness of the observed community. Only 2 RFLP pro- 
files (RHIZO-1 and RHIZO-2) were found in both samples. The 
whole root communities were more uniform with as many as 
five shared RFLP profiles. In Rl (25 observed individuals) and 
R2 (22 observed individuals), RHIZO-1 was the most abun- 
dant RFLP profile (48 and 50%, respectively). Another 4 RFLP 
profiles were shared, ACAU-1 (16 and 18.2%), ARCH-1 (8 and 
18.2%), ARCH-2 (8 and 9.1%), and RHIZO-2 (12 and 4.5%). 
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FIGURE 3 | Typical AMF structures, stained with 0.1% cotton blue, 
encountered in C. japonica 'Nuccio's Pearl' roots. Intraradical hyphae (IH), 
along with intraradical vesicles (V) that were recognized because of the 



presence of lipid droplets (LD), pervaded the root apparatuses. Intracellular 
arbuscules (A) colonized large patches of the root cortex, as expected, 
without entering the central cylinder (CC). The bars stand for 30 |im. 
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FIGURE 4 | Percentage of aspecific amplification detected in the whole 
root samples (R1 and R2) and from 100 (A100), 300 (A300), and 500 
(A500) arbusculated cell aliquots, respectively. 
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FIGURE 5 | RFLP profile percentage distribution for the AMF 
communities in the whole root samples (R1 and R2) and in the 300 
(A300) and 500 (A500) arbusculated cell aliquots. 



One RFLP profile, DIVE-2, was only found in Rl (8%). ESC 
was 100% for Rl and 95.5% for R2. Bias corrected CHAOl 
gave 6.0 ± 0.0 and 5.0 ± 0.3 expected RFLP profiles for Rl and 
R2, respectively. Only for sample R2, was the asymptote not 
reached by the related rarefaction curve (Figure 6). Overall, the 
cumulative AMF communities from the whole root (WR) and 
LMD AC shared 4 RFLP profiles, RHIZO-1, ARCH-1, ARCH- 
2, and RHIZO-2 (Figure 7). In this case, the rarefaction curves 
related to both cumulative communities reached the asymptote 
(Figure 6). 

DISCUSSION 

To the authors' knowledge, this experiment has been the first 
attempt to characterize an arbuscule forming AMF community 
through the use of LMD. Previous studies on AM symbiosis using 



LMD were aimed at investigating gene expression on AMF iso- 
lates that were inoculated singularly in microcosm conditions 
(Balestrini et al., 2007; Guether et al, 2009; Giovannetti et al, 
2012). In the present study, a mesocosm approach has been 
adopted. The used inoculum originated from natural soils and 
it is therefore likely that the AMF species that were present were 
propagated along with other rhizosphere biota. 

Prior to DNA extraction, three aliquots, containing 100, 300, 
and 500 AC (A100, A300, and A500, respectively), were pre- 
pared by dissecting them from root sections from at least 10 
different root fragments. The nested PCR approach, using the 
primer pair AML1-AML2, gave a positive band of averagely 
800 bp for all three samples, as expected. After cloning and RFLP 
typing, a sequencing step made it possible to characterize the 
taxonomical origin of the DNA amplicons. All the sequences 
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FIGURE 6 | Rarefaction curves indicating RFLP profile accumulation 
related to the number of screened clones for the LMD samples (A300, 
A500), whole root samples (R1, and R2) and respective cumulative 



samples (AC and WR). The dotted lines indicate the upper and lower 
95%-confidence intervals. The red lines indicate that the curve reached the 
asymptote. 




FIGURE 7 | RFLP profile percentage distribution for the whole root 
cumulative community (R1 + R2 = WR) and the arbusculated cell 
cumulative community (A300 + A500 = AR). 



retrieved from sample A100 were the result of aspecific ampli- 
fications. This was possibly due to the fact that this primer 
pair (Lee et al., 2008) as well as other primer pairs in gen- 
eral, although said to be Glomeromycota-specific, can sometimes 
lead to aspecific amplification when the starting template resides 
in a complex matrix and is overwhelmed by non-target DNA 
(Douhan et al, 2005; Alguacil et al, 2008; Toljander et al, 2008; 
Borriello et al., 2012). Samples A300 and A500 gave decreasing 
rates of non-target amplification, while aspecific amplification 
was absent or negligible for the whole root samples (Rl and 
R2). Around 50% of the screened clones in sample A300 were 
shown to be carrying a non-target fragment and this led to 
a decreased number of total observed AMF individuals (13). 
This factor was clearly responsible for the lower sample cover- 
age achieved for this sample (ESC = 91.7%, with the rarefaction 



curve not reaching the asymptote, Figure 6). The reciprocal dis- 
crepancies observed in terms of RFLP profile presence/abundance 
(RHIZO-1, ARCH-1, ARCH-2, and GLO-6) between A300 and 
A500 suggest that these two samples singly cannot be regarded 
as 100% representative of the actual arbuscule forming commu- 
nity. This could be due to a lack of robustness in the sampling 
effort, although the coverage and species richness estimators 
denoted a remarkable representativeness of the sample com- 
munities. It is possible that the complexity of the mesocosm 
conditions may have triggered the development of a heteroge- 
neous rhizosphere biota. This could have played a critical role 
in diminishing or making AMF detection impossible in smaller 
cell aliquots. The spectrophotometer measurements of the DNA 
extraction yields from the different cell aliquots showed com- 
parable values and did not differ remarkably from the ones 
performed on the whole root extractions. However, these val- 
ues refer to total DNA and may not reflect the actual presence 
of target AMF DNA in the sample. Moreover, AC were col- 
lected on a morphological basis and, to avoid an impact on DNA 
preservation, no staining was performed to aid their discrim- 
ination from other non-AMF structures. Thus, further studies 
will require a greater effort in terms of number of AC collected 
from roots sampled from mesocosm or natural conditions. As 
an alternative, new and more specific primer pairs could be 
designed and their use consolidated, in order not to lose AMF- 
specificity due to the overwhelming pressure exerted by non- 
target organisms that are present in complex matrices. The use of 
DNA conservative fungal staining techniques (Pitet et al., 2009) 
should also be considered, in order to assist arbusculated cell 
individuation. 

In general, the interaction between the involved symbiotic 
partners was significant as nearly every mycorrhizal root part 
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checked presented a profuse arbuscule production. This, in 
combination with a high frequency and intensity of colonization, 
made it easy to collect suitable root fragments for LMD use, 
as the technique relies mainly on the degree of availability and 
morphological detectability of intracellular arbuscules. Overall, 
out of the 27 taxa inoculated in the substrate (see Materials and 
Methods), only 5 were recovered from microdissected AC. This is 
an important result, as it could suggest that the related 5 taxa are 
capable of undertaking a close partnership with the plant species 
under study and which actively express their functional traits in 
terms of benefits for plant fitness. A large number of ultraradi- 
cal vesicles were also found all along the root system (Figure 3). 
Vesicles can carry a comparable number of nuclei to that of 
the spores (Camper et al., 2008) and can therefore be detected 
through PCR with specific primers. This may lead to an overde- 
tection of the fungal isolates that produce abundant intraradical 
vesicles, when the AMF community composition of the whole 
root is being investigated. This could be the case of the fun- 
gal isolate affiliated to RFLP profile RHIZO-2, which clustered, 
according to the phylogenetic tree, with a sequence belonging to 
Rhizophagus fasciculatum (Figure 2) and was the most abundant 
type found in the roots (Figure 7). This AMF species is char- 
acterized by profuse vesicle formation (Gerdemann and Trappe, 
1974) and the greater presence of gDNA belonging to this taxon 
in the root samples might therefore have hindered or precluded 
the molecular recovery of other coexisting taxa. This is the case of 
ARCH-1 and ARCH-2, which were less abundant in the whole 
root community with respect to the community that featured 
AC. These RFLP profiles grouped with sequences affiliated to the 
Archaeosporales order, a taxonomical rank that includes species 
that form rare or no vesicles (Morton and Redecker, 2001; Walker 
et al., 2007). The same considerations can be made for GLO- 
6, which was not found in the whole root samples but was, 
albeit scarcely, detectable in the community from the LMD AC. 
Possible reasons for this could be that the small patches infected 
by this fungus along the roots were not spotted after whole root 
DNA extraction, due to the overwhelming presence of DNA of 
other more abundant taxa, but were instead detected through 
LMD. The explanation could be that this technique, being based 
only on the collection of intracellular arbuscules, can eliminate 
the co-extraction of DNA from other fungal structures (vesicles, 



spores, hyphae) that exert a dilution effect on the actual arbuscule 
forming community. Therefore, it is likely that, thanks to LMD, 
candidate beneficial AMFs that colonize small patches within the 
root, but which play a pivotal role in exploring the surrounding 
soil (Parniske, 2008), can be spotted. Interestingly, two RFLP pro- 
files, i.e., ACAU- 1 and DIVE-2, were highlighted in the whole root 
samples but not in the LMD samples, indicating that they could 
be located within the root but might not produce intracellular 
arbuscules. This finding would seem to support the hypothesis 
that the production of arbuscules by different AMF taxa could 
be differentially stimulated according to their momentary state 
of efficiency. Conceptually, it is feasible that ACAU-1 and DIVE- 
2 could be AMF species that were momentarily expressing less 
beneficial attributes for the plant at the moment of sampling. 
Interestingly, recent findings have shown how AM symbiosis, 
depending on the partners involved, can range from coopera- 
tive to less-cooperative (Verbruggen and Kiers, 2010; Kiers et al., 
2011). However, to provide supporting experimental evidence, 
further analyses are necessary to investigate the functional traits 
involved in AM symbiosis. 

In conclusion, the LMD technique can be considered a power- 
ful tool for the obtainment of a more qualitative representation of 
the arbuscule forming and, therefore, symbiotically active AMF 
community. Future studies on LMD-mediated identification of 
active AMF should focus on the field validation of this method 
for AMF traceability studies and could be finalized on the formu- 
lation of selected inocula for field and greenhouse applications. 
Further efforts should be devoted to LMD-mediated functional 
biodiversity studies, in which community composition analyses 
based on functional genes, e.g., the phosphate transporter gene, 
will be able to shed light on the beneficial traits involved in AM 
symbiosis. 

ACKNOWLEDGMENTS 

The authors wish to thank Fiori Tipici del Lago Maggiore e del 
Biellese S.c.r.l. and Tecnoverde di Cesa S.p.a. (Verbania, Italy) 
for the technical support, as well as R Bizioli, G Morandi, 
and MT Delia Beffa for their assistance in the experimental 
work; ECOFLOR project (Regione Piemonte) and CNR (Premiale 
2011 project) for funding. AB's fellowship was funded by the 
PRO-LACTE project (Regione Piemonte). 



REFERENCES 

Alguacil, M. M., Lumini, E., Roldan, 
A., Salinas-Garcia, J. R., Bonfante, 
R, and Bianciotto, V. (2008). The 
impact of tillage practices on arbus- 
cular mycorrhizal fungal diversity 
in subtropical crops. Ecol. Appl 18, 
527-536. 

Balestrini, R., and Bonfante, P. 
(2008}. Laser microdissection 
(LM): applications to plant 
materials. Plant Biosyst. 142, 
331-336. 

Balestrini, R., Goimez-Ariza, J., 
Lanfranco, L., and Bonfante, 
P. (2007). Laser microdissec- 
tion reveals that transcripts for 



five plant and one fungal phos- 
phate transporter genes are 
contemporaneously present in 
arbusculated cells. Mol. Plant 
Microbe Interact. 20, 1055-1062. 

Balestrini, R., Magurno, R, Walker, 
C, Lumini, E., and Bianciotto, 
V. (2010). Cohorts of arbuscular 
mycorrhizal fungi (AMF) in Vitis 
vinifera, a typical Mediterranean 
fruit crop. Environ. Microbiol. Rep. 
2, 594-604. 

Bordello, R., Lumini, E., Girlanda, 
M., Bonfante, P., and Bianciotto, V. 
(2012). Effects of different manage- 
ment practices on arbuscular myc- 
orrhizal fungal diversity in maize 



fields by a molecular approach. Biol. 
Fert. Soils 48, 911-922. 

Chao, A., and Shen, T. J. (2010). 
Program SPADE (Species Prediction 
and Diversity Estimation). 
Available online at: http://chao. 
stat.nthu.edu.tw 

Cole, J. R., Chai, B., Marsh, T. L., 
Farris, R. J., Wang, Q., Kulam, 
S. A., et al. (2003). The ribo- 
somal database project (RDP-II): 
previewing a new autoaligner that 
allows regular updates and the new 
prokaryotic taxonomy. Nucl. Acids 
Res. 31,442-443. 

Davison, J., Opik, M., Zobel, M., Vasar, 
M., Metsis, M., and Moora, M. 



(2012). Communities of arbus- 
cular mycorrhizal fungi detected 
in forest soil are spatially het- 
erogeneous but do not vary 
throughout the growing season. 
PLoS ONE 7:e41938. doi: 10.1371/ 
journal.pone.0041938 
Day, R. C, Grossniklaus, U., and 
MacKnight, R. C. (2005). Be 
more specific! Laser-assisted 
microdissection of plant 
cells. Trends Plant Set. 10, 
397-405. 

Douhan, G. W., Petersen, C, Bledsoe, 
C. S., and Rizzo, D. M. (2005). 
Contrasting root associated fungi of 
three common oak-woodland plant 



Frontiers in Plant Science | Plant-Microbe Interaction 



May 2013 | Volume 4 | Article 135 | 8 



Berruti et al. 



Laser microdissection for AMF detection 



species based on molecular iden- 
tification: host specificity or non- 
specific amplification? Mycorrhiza 
15, 365-372. 
Dumbrell, A. J., Ashton, P. D., 
Aziz, N., Feng, G., Nelson, 
M., Dytham, C, et al. (2011). 
Distinct seasonal assemblages 
of arbuscular mycorrhizal fungi 
revealed by massively parallel 
pyrosequencing. New Phytol. 190, 
794-804. 

Felsenstein, J. (1985). Confidence lim- 
its on phylogenies: an approach 
using the bootstrap. Evolution 39, 
783-791. 

Gamper, H. A., Young, J. P. W., Jones, 
D. L., and Hodge, A. (2008). Real- 
time PCR and microscopy: are the 
two methods measuring the same 
unit of arbuscular mycorrhizal fun- 
gal abundance? Fungal Genet. Biol 
45, 581-596. 

Gaude, N., Bortfeld, S., Duensing, 
N., Lohse, M., and Krajinski, F. 
(2011). Arbuscule-containing and 
non-colonized cortical cells of myc- 
orrhizal roots undergo a massive 
and specific reprogramming 
during arbuscular mycor- 
rhizal development. Plant J. 69, 
510-528. 

Gerdemann, J. W., and Trappe, J. M. 
(1974). The endogonaceae in the 
Pacific Northwest. Mycol. Mem. 5, 
1-76. 

Giovannetti, M., Balestrini, R., Volpe, 
V. , Guether, M . , Straub, D. , 
Costa, A., et al. (2012). Two 
putative- aquaporin genes are 
differentially expressed during 
arbuscular mycorrhizal sym- 
biosis in Lotus japonicus. BMC 
Plant Biol. 12:186. doi: 10.1186/ 
1471-2229-12-186 

Gomez, S. K., Javot, H., 
Deewatthanawong, P., Torres- 
Jerez, I., Tang, Y., Blancaflor, E. 
B., et al. (2009). Medicago trun- 
catula and Glomus intraradices 
gene expression in cortical cells 
harboring arbuscules in the arbus- 
cular mycorrhizal symbiosis. BMC 
Plant Biol. 9:10. doi: 10.1186/ 
1471-2229-9-10 

Gupta, R. K., and Sharma, C. (2010). 
Diversity of arbuscular mycor- 
rhizal fungi in Camellia sinensis 
in Uttarakhand State, India. Afr. J. 
Biotechnol. 9,5313-5319. 

Guether, M., Balestrini, R., Hannah, 
M., He, J., Udvardi, M., and 
Bonfante, P. (2009). Genome-wide 
reprogramming of regulatory 
networks, transport, cell wall 
and membrane biogenesis during 
arbuscular mycorrhizal symbiosis 
in Lotus japonicus. New Phytol. 182, 
200-212. 



Hammer, 0., Harper, D. A. T, 
and Ryan, P. D. (2001). PAST: 
paleontological statistics software 
package for education and data 
analysis. Palaeontol. Electron. 4, 1-9. 

Hogekamp, C, Arndt, D., Pereira, 
P., Becker, J. D., Hohnjec, N., 
and Kiister, H. (2011). Laser- 
microdissection unravels cell-type 
specific transcription in arbuscu- 
lar mycorrhizal roots, including 
CAAT-box TF gene expression 
correlating with fungal contact 
and spread. Plant Physiol. 157, 
2023-2043. 

Katoh, K., and Toh, H. (2008). Recent 
developments in the MAFFT 
multiple sequence alignment 
program. Brief. Bioinform. 9, 
286-298. 

Kiers, E. T., Duhamel, M., Beesetty, 
Y., Mensah, J. A., Franken, O., 
Verbruggen, E., et al. (2011). 
Reciprocal rewards stabilize 
cooperation in the mycor- 
rhizal symbiosis. Science 333, 
880-882. 

Kimura, M. (1980). A simple method 
for estimating evolutionary rate 
of base substitutions through 
comparative studies of nucleotide 
sequences. /. Mol. Evol. 16, 111-120. 

Kriiger, M., Kriiger, C, Walker, C., 
Stockinger, H., and Schii61er, 
A. (2012). Phylogenetic refer- 
ence data for systematics and 
phylotaxonomy of arbuscular 
mycorrhizal fungi from phylum 
to species-level. New Phytol. 193, 
970-984. 

Lee, J., Lee, S., and Young, J. P. W. 
(2008). Improved PCR primers 
for the detection and identifica- 
tion of arbuscular mycorrhizal 
fungi. FEMS Microbiol. Ecol. 65, 
339-349. 

Liu, R. J., and Wang, F. Y. (2003). 
Selection of appropriate host plants 
used in trap culture of arbuscular 
mycorrhizal fungi. New Phytol. 115, 
495-501. 

Lumini, E., Orgiazzi, A., Bordello, 
R., Bonfante, P., and Bianciotto, 
V. (2010). Disclosing arbuscular 
mycorrhizal fungal biodiversity 
in soil through a land-use gra- 
dient using a pyrosequencing 
approach. Environ. Microbiol. 12, 
2165-2179. 

Lumini, E., Vallino, M., Alguacil, M. 
M., Romani, M., and Bianciotto, 
V. (2011). Different farming and 
water regimes in Italian rice fields 
affect arbuscular mycorrhizal fun- 
gal soil communities. Ecol. Appl. 21, 
1696-1707. 

Morton, J. B., and Redecker, D. 
(2001). Two new families of 
Glomales, Archaeosporaceae and 



Paraglomaceae, with two new gen- 
era Archaeospora and Paraglomus, 
based on concordant molecular 
and morphological characters. 
Mycologia 93, 181-195. 
Nelson, T., Tausta, S. L., Gandotra, 
N., and Liu, T. (2006). Laser 
microdissection of plant tis- 
sue: what you see is what you 
get. Ann. Rev. Plant Biol. 57, 
181-201. 

Oehl, E, Sieverding, E., Mader, P., 
Dubois, D., Ineichen, K., Boiler, T., 
et al. (2004). Impact of long-term 
conventional and organic farm- 
ing on the diversity or arbuscular 
mycorrhizal fungi. Oecologia 138, 
574-583. 

Oehl, E, Sieverding, E., Palenzuela, 
J., Ineichen, K., and Alves da 
Silva, G. (2011). Advances in 
Glomeromycota taxonomy and 
classification. IMA Fungus 2, 
191-199. 

Opik, M., Metsis, M., Daniell, T J., 
Zobel, M., and Moora, M. (2009). 
Large-scale parallel 454 sequenc- 
ing reveals host ecological group 
specificity of arbuscular mycor- 
rhizal fungi in a boreonemoral 
forest. New Phytol. 184, 424-^37. 

Opik, M., Zobel, M., Cantero, J. J., 
Davison, J., Facelli, J. M., Hiiesalu, 
L, et al. (2013). Global sam- 
pling of plant roots expands the 
described molecular diversity of 
arbuscular mycorrhizal fungi. 
Mycorrhiza. doi: 10.1007/s00572- 
013-0482-2. [Epub ahead of print]. 

Pandey, A., and Palni, L. M. S. 
(1996). The rhizosphere effect 
of tea on soil microbes in a 
Himalayan monsoonal loca- 
tion. Biol. Fertil. Soils 21, 
131-137. 

Parniske, M. (2008). Arbuscular myc- 
orrhiza: the mother of plant root 
endosymbioses. Nat. Rev. Microbiol. 
6, 763-775. 

Pellegrino, E., Turrini, A., Gamper, 
H. A., Cafa, G., Bonari, E., Young, 
J. P., et al. (2012). Establishment, 
persistence and effectiveness of 
arbuscular mycorrhizal fungal 
inoculants in the field revealed 
using molecular genetic tracing and 
measurement of yield components. 
New Phytol. 194, 810-822. 

Perez-Tienda, J., Testillano, P. 
S., Balestrini, R., Fiorilli, V., 
Azcon-Aguilar, C., and Ferrol, 
N. (2011). GintAMT2, a new 
member of the ammonium trans- 
porter family in the arbuscular 
mycorrhizal fungus Glomus 
intraradices. Fungal Genet. Biol. 48, 
1044-1055. 

Pitet, M., Camprubi, A., Calvet, C., 
and Estaun, V. (2009). A modified 



staining technique for arbuscu- 
lar mycorrhiza compatible with 
molecular probes. Mycorrhiza 19, 
125-131. 

Porter, W. M. (1979). The "most prob- 
able number" method for enumer- 
ating infective propagules of vesic- 
ular arbuscular mycorrhizal fungi 
in soils. Aust. J. Soil Res. 17, 
515-519. 

Saitou, N., and Nei, M. (1987). 
The neighbor-joining method: a 
new method for reconstructing 
phylogenetic trees. Mol. Biol. Evol. 

4, 406-425. 

Smith, S. E., and Read, D. J. (2008). 

Mycorrhizal Symbiosis. 3rd Edn. 

Amsterdam: Academic Press. 
Tamura, K., Peterson, D., Peterson, N., 

Stecher, G., Nei, M., and Kumar, 

5. (2011). MEGA5: molecular 
evolutionary genetics analysis using 
maximum likelihood, evolutionary 
distance, and maximum parsi- 
mony methods. Mol. Biol. Evol. 28, 
2731-2739. 

Tisserant, E., Kohler, A., Dozolme- 
Seddas, P., Balestrini, R., 
Benabdellah, K., Coiard, A., 
et al. (2012). The transcriptome 
of the arbuscular mycorrhizal 
fungus Glomus intraradices (DAOM 
197198) reveals functional trade- 
offs in an obligate symbiont. New 
Phytol. 193, 755-769. 

Toljander, J. E, Santos- Gonzalez, J. C, 
Tehler, A., and Finlay, R. D. (2008). 
Community analysis of arbuscu- 
lar mycorrhizal fungi and bac- 
teria in the maize mycorrhizo- 
sphere in a long-term fertilization 
trial. FEMS Microbiol. Ecol. 65, 
323-338. 

Trouvelot, A., Kough, J., and 
Gianinazzi-Pearson, V. (1986). 
"Evaluation of VA infection levels in 
root systems. Research for estima- 
tion methods having a functional 
significance," in Physiological and 
Genetical Aspects of Mycorrhizae 7 
eds V. Gianinazzi-Pearson, S. 
Gianinazzi (France: INRA Press), 
217-221. 

Verbruggen, E., and Kiers, E. T. (2010). 
Evolutionary ecology of mycorrhizal 
functional diversity in agricultural 
systems. Evol. Appl. 3, 547-560. 

Walker, C, Vestberg, M., Demircik, 
E, Stockinger, H., Saito, M., 
Sawaki, H., et al. (2007). Molecular 
phylogeny and new taxa in the 
Archaeosporales (Glomeromycota): 
Ambispora fennica gen. sp nov., 
Ambisporaceae fam. nov., and 
emendation of Archaeospora and 
Archaeosporaceae. Mycol. Res. Ill, 
137-153. 

White, T. J., Bruns, T, Lee, S., and 
Taylor, J. (1990). "Amplification 



www.f ro ntiersin.org 



May 2013 | Volume 4 | Article 135 | 9 



Berruti et al. 



Laser microdissection for AMF detection 



and direct sequencing of fungal 
ribosomal RNA genes for phy- 
logenetics" in PCR Protocol: A 
Guide to Methods and Applications, 
eds M. A. Innis, D. H. Gelfand, 
J. S. Sninsky, and T. J. White 
(San Diego, CA: Academic Press), 
315-322. 

Zhang, Z., Schwartz, S., Wagner, L., 
and Miller, W (2000). A greedy 
algorithm for aligning DNA 
sequences. /. Comput. Biol. 7, 
203-214. 



Zheng, J., Peneva, V., and Brown, D. J. 
F, (2000). Longidorus camelliae n. 
sp. (Nematoda: Longidoridae) asso- 
ciated with ornamental cultivars 
of Camellia japonica L. growing 
in a nursery at Fuyang, Zhejiang 
Province, eastern China. Syst. 
Parasitol. 47, 119-125. 

Conflict of Interest Statement: The 

authors declare that the research 
was conducted in the absence of any 
commercial or financial relationships 



that could be construed as a potential 
conflict of interest. 

Received: 15 February 2013; accepted: 22 
April 2013; published online: 09 May 
2013. 

Citation: Berruti A, Bordello R, Lumini 
E, Scariot V, Bianciotto V and Balestrini 
R (2013) Application of laser microdis- 
section to identify the mycorrhizal fungi 
that establish arbuscules inside root cells. 
Front. Plant Sci. 4:135. dot: 10.3389/fpls. 
2013.00135 



This article was submitted to Frontiers in 
Plant-Microbe Interaction, a specialty of 
Frontiers in Plant Science. 
Copyright © 2013 Berruti, Borriello, 
Lumini, Scariot, Bianciotto and 
Balestrini. This is an open-access 
article distributed under the terms of the 
Creative Commons Attribution License, 
which permits use, distribution and 
reproduction in other forums, provided 
the original authors and source are cred- 
ited and subject to any copyright notices 
concerning any third-party graphics etc. 



Frontiers in Plant Science | Plant-Microbe Interaction 



May 2013 | Volume 4 | Article 135 | 10 



